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We describe a semi-empirical atomic basis Extended Hiickel Theoretical (EHT) technique that 
can be used to calculate bulk bandstructure, surface density of states, electronic transmission and 
interfacial chemistry of various materials within the same computational platform. We apply this 
method to study multiple technologically important systems, starting with carbon-nanotubes (CNT) 
and their interfaces in this paper, and silicon-based heterostructures in our follow-up paper. We 
find that when it comes to quantum transport through interesting, complex heterostructures, the 
Hiickel bandstructure offers a fair and practical compromise between orthogonal tight-binding the- 
ories (0TB) with limited transferability between environments under large distortion, and density 
functional theories (DFT) that are computationally quite expensive for the same purpose. 

PACS numbers: 61.46.Fg, 73.43. Cd, 73.63.Fg 
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I. INTRODUCTION 



Quantitative electronic structure theories are essential 
to the understanding and designing of novel materials 
and devices. It is now generally accepted that transport 
properties of nanoscale devices depend on both the in- 
trinsic electronic structure of the active channel, as well 
as its interfacial properties with contacts and other scat- 
tering centers. A particular challenge in this respect is to 
incorporate both long and short range correlations within 
the same framework, such as the bulk bandstructure of 
periodic solids as well as the local chemical properties 
of clusters, surfaces, and interfaces. For instance, simu- 
lating scanning tunneling spectra (STS) of molecules on 
silicon substrates requires an accurate description of the 
silicon bulk bandstructure that quantitatively captures 
not just the bandgap responsible for the onset of nega- 
tive differential resistance^i but also the multiple effective 
masses which determine the contact density of states and 
injection velocities, and the strain parametrizations that 
capture atomic reconstruction and relaxation near the 
surface and their bonding with molecular componentsj^ 
In addition, one needs to describe the electrostatics re- 
sponsible for band-bending in the silicon depletion layer, 
the molecular transport levels and their transmission un- 
der bias, and finally the electronic properties of the scan- 
ning tip and the intervening vacuum layer, all within the 
same formalism. It is no surprise therefore that under 
these circumstances, standard electronic structure tech- 
niques developed by quantum chemists for simulating 
molecules are usually incompatible with those developed 
by solid-state physicists for bulk bandstructure, making 
it important to develop a common formalism that ad- 
dresses both domains of interest and also maintains a 
good compromise between computational accuracy and 
practicality. 

While sophisticated methods exist for equilibrium ge- 



ometry and bandstructure, comparable success has yet 
to be achieved for transport problems, partly because of 
the lack of universally accepted experimental standards, 
but mainly because quantum transport inherently in- 
volves solving a complicated nonequilibrium open bound- 
ary problem for which electronic structure theories are 
not benchmarked. A proper quantitative understanding 
of correlation effects in transport is still evolving and it is 
not yet clear if mean-field approaches that work at equi- 
librium are at all capable of handling the profusion of 
electronic excitations that often dominate in nanoscale 
conduction^i4 Even aside from such correlation issues, 
one needs to worry about heterointerfaces since current 
flow involves charge transport across two intrinsically dif- 
ferent material systems - a multimoded contact consist- 
ing of a highly conductive material externally maintained 
at thermal equilibrium, and a sparsely moded device re- 
gion that is readily driven away from equilibrium and 
acts as the active transport channel. 

In this paper, we employ a semi-empirical approach to 
electronic structure that can be adopted for electronic 
conduction through complex hybrid systems by combin- 
ing it with the Nonequilibrium Green's function (NEGF) 
technique for quantum transport. Our theoretical pa- 
rameters are tailored to salient features of the bulk band- 
structure, while the employment of non-orthogonal basis 
sets resembling linear combinations of underlying atomic 
orbitals seems to make them fairly transferrable to sur- 
faces as well, as observed in the pastil In addition, 
the presence of explicit basis sets opens up the possi- 
bility of 'stitching' together disparate regions^ by match- 
ing the interfacial Green's function, which is the only 
quantity through which the diverse regions communicate 
with each other quantum mechanically. The modularity 
of our approach also allows us to conveniently separate 
the problems of determining the optimized interfacial ge- 
ometry and the interfacial transmission (we are ignoring 
current-induced forces), the former depending on the to- 
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tal energy of the system while the latter depends only 
on a few relevant single-electron levels near the Fermi 
energy. In other words, given a particular atomistic con- 
figuration of the contact- channel- contact hetero structure, 
we seek to determine its transport properties by coupling 
our electronic structure approach with quantum transport 
using NEGF. 

The outline of the paper is as follows: section II ex- 
plains the strengths of EHT over other traditional band- 
structure methods. In section III we briefly summarize 
the NEGF aproach used to calculate density of states 
and transmission of CNTs; we then specify the model 
Hamiltonian and describe the details of EHT used to de- 
termine the bandstructure. The numerically calculated 
bandstructure data for nanotubes are then compared in 
section IV with experimental scanning tunneling spec- 
troscopy experiments along with other theoretical ap- 
proaches. In section V we investigate the changes in the 
dispersion of a semi-conducting CNT under large lateral 
deformations as well as with a CO molecular attachment 
to its surface that allows it to function as a molecular 
sensor. We summarize our work and discuss future ex- 
tensions in the last section. 



II. WHY THIS PARTICULAR METHOD? 

A particular trade-off in any bandstructure theory is 
between flexibility and rigor. While empirical, orthogo- 
nal tight-binding (0TB) methods are quick and practi- 
cal, they are bcnchmarked for specific geometries and arc 
usually not very transferable to other environments in- 
volving significant structural deformations beyond a few 
percent. Tight binding basis-sets are commonly assumed 
to be both orthogonal and short-ranged while atomic 
wavefunctions are not, meaning that 0TB basis sets do 
not resemble eigenstates of an atomic Hamiltonian. Ef- 
forts at improving tight-binding theories involve going 
beyond nearest neighbor techniques, using higher virtual 
orbital bases for increased completenessj^iSiiS and em- 
ploying power laws for parameter transfer under small 
(~ 1 — 2%) straiuiii Nevertheless these models are likely 
to miss important chemical details involving properties 
varying on an interatomic length scale, in particular near 
deformed surfaces and interfaces where a drastic recon- 
struction of the atomic structure is expectedi^ 

At the other end of the spectrum are accurate, but 
computationally expensive first-principles techniques de- 
veloped by quantum chemists and solid state physicists, 
such as Configuration Interaction (CI) and Density Func- 
tional Theories (DFT) in various atomic or plane wave 
basis sets or combinations thereof. Structural deforma- 
tions are naturally captured by such total energy cal- 
culations by solving a one electron Schrodinger equa- 
tion in a suitable self-consistent potential approximat- 
ing the electron-electron interaction>i2iiii4 Such codes 
are typically based on rigorous variational theorems and 
are quantitatively quite accurate, at least for equilib- 



rium properties. Their extension and practical im- 
plementation to transport beyond the linear response 
regime is continuously evolving ^^iiSiii and a topic of cur- 
rent researchii^ Conceptually, it is not clear if any self- 
consistent potential approach can quantitatively describe 
the rich spectrum of many-body transitions that arc often 
experimentally accessed in strongly correlated transport 
in weakly coupled systems^*^ 

We aim for a practical compromise between these two 
limits by using a semi-empirical technique motivated by 
Extcnded-Hiickel calculations popular in the chemistry 
community. Such theories, widely used in the past to de- 
scribe the equilibrium properties of isolated molecules^ 
have recently been applied to molecular conductionSS and 
also extended to solids using transferable atomic-orbital 
basis sets (AO) for calculating the electronic structure 
of various compounds benchmarked with detailed DFT 
calculations within the local density (LDA) or general- 
ized gradient (GGA) approximations^ Given a geom- 
etry, one uses the explicit EHT basis functions to cal- 
culate a non-orthogonal overlap matrix S", which along 
with separately fitted onsite Hamiltonian matrix ele- 
ments yields the corresponding off-diagonal hopping ele- 
ments of the Hamiltonian. Within the standard Hiickcl 
prescription, structural changes are simply accounted for 
by re-calculating the overlap and hopping elements, but 
leaving the basis sets and onsite elements unchanged. 

In the following, we apply this EHT parametriza- 
tion schemeSi by benchmarking it to a two-dimensional 
graphene sheet^SiS^ and extending it to obtain the band- 
structure, density of states and electronic transmission 
of carbon nanotubes (CNT) of varying chiralities. We 
show that the same bulk-optimized EHT-paramcters (on- 
site energies and AO-basis functions) are transferable 
to small diameter CNT bandstructures, capturing even 
curvature-induced bandgap effects for larger than 1 — 3% 
tube deformation, in quantitative agreement with STS 
data. Furthermore, surface chemical effects are examined 
through the study of nanotube based carbon-monoxide 
sensors whose alteration of electronic structure upon 
molecular adsorption compares favorably with ab-initio 
calculations of da Silva et al.m^ In our follow-up pa- 
per, we will demonstrate a similar transferability be- 
tween bulk silicon, various silicon surfaces, apply the 
EHT-methodology to unreconstructed silicon nanowires. 
Taken together, the wide variety of these examples illus- 
trates the range of transferability of EHT parameters, 
which we believe makes Extended Hiickcl Theory a use- 
ful practical tool for electronic structure and quantum 
transport. 



III. APPROACH 

Simulating conduction through a heterostructure in- 
volves combining suitable bandstructures for the channel 
and contact materials with self-consistent electrostatics 
and quantum transport. While the formulation of cor- 
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related transport is itself an active area of research, our 
aim here is to develop a minimal model that would cap- 
ture quantum chemistry, surface physics, bandstructure 
and electrostatic effects that are all crucial for different 
prominent aspects of nanoscale conduction. The ingredi- 
ents needed for a proper simulation are the Hamiltonian 
and potential matrices describing the device bandstruc- 
ture and electrostatics, and the contact self-energies that 
effectively open up the system and allow us to add or re- 
move charge under nonequilibrium conditions. The EHT 
prescription gives us a practical way to calculate these 
ingredient matrices for a given atomistic structure, and 
then connect them with a non-equilibrium Green's func- 
tion (NEGF) formulation of quantum transport, which 
we briefly summarize below^S*^ The retarded Green 
function of the device is given by 

G^[{E + ii^)S -h-j:i- S,]"' , (1) 

where S and H describe the overlap and Hamiltonian 
matrices of the device unit cell calculated according to 
the Hiickel prescription Eq.©. The matrix elements S/ 
are self-energies that provide open boundary conditions 
to the device with the left and right semi-infinite con- 
tacts. The self-energy S = rgr'^ incorporates the cou- 
pling matrix t describing the contact-device bonding, 
while g is the surface Green function of the left/right con- 
tact calculated by means of a recursion technique 
In a non-orthogonal tight-binding scheme the density 
of states (DOS) is given by D{E) = ^Tr {AS) where 
A = i(G — G^) denotes the spectral function. Finally, 
in the phase-coherent limit the zero-bias transmission 
through the unit cell reads T[E) = Tr [TiGTrC^] , where 
Ti^r = i {^i,r — Sz,r) are the broadening matrices which 
specify the time an electron resides within the device. 
In this paper, we will study infinite nanotubes so that 
the active device is just one CNT unit cell and the left 
and right contacts extend that cell to infinity in either 
direction. 

The bandstructure of a nanotube with chirality 
(n, m) is calculated employing the non-orthogonal Slater- 
Koster scheme and solving for the generalized eigenvalue 
problemSS 

H(k)^,(k) =^,(k)S(k)A,(k) , (2) 

where ^^(k) denotes the eigenvector of the i*^ subband, 
and k is a Bloch wavevector within the first Brillouin 
zone. The size of the overlap and Hamiltonian matrices 
are determined by the chosen basis set, i.e., the number 
of atoms within the unit cell multiplied by the number 
of orbitals per atom. In our case, using four sp-orbitals 
per atom, the size of these matrices is 80 x 80 for an 
armchair (5, 5) tube as sketched in Fig. ^ The overlap- 
and Hamilton matrices S(k) and H(k) representing the 
structure in reciprocal space are calculated by 

H,,,(k) = e'''-(^--^.''"') H,o.,'™' (3) 

j' ,m' 



S,,,(k) = e^k.(R.o-R,,„0 ^ (4) 

j' .m' 

where i,j label atoms within the unit cell, and m' is the 
unit cell index. The summation indices in H^j- (k). Si j- (k) 
run over all atoms j' in unit cell m! which are equivalent 
to atom j in the reference unit cell m = 0. The real-space 
matrix elements Hio.j'm' , SiOj'm', between an atom i 
in the reference unit cell and atom j' in cell m' are cal- 
culated by means of the Extended Hiickel prescription. 




FIG. 1: Sketch of a (5, 5) armchair tube with three unit cells. 
The dashed rectangle marks the center unit cell containing 20 
carbon atoms. For the calculation of the E — k dispersion we 
use the 1st nearest-neighbor unit cells to the left and right, 
respectively. 

Perhaps the most prominent difference between the 
empirical tight-binding and EHT principle is the use of 
an irreducible representation of the tight-binding Hamil- 
tonian which is directly fitted to available bandstruc- 
ture data without employing explicit basis sets. In EHT 
however one works directly with the orbital basis func- 
tions out of which the Hamiltonian elements are con- 
structed using the Hiickel principle. The diagonal ele- 
ments are benchmarked with experimental values of elec- 
tronic 'hardness', i.e. the difference between ionization 
potential and electron affinityiiS The off-diagonal matrix 
elements are then determined directly from the prescrip- 
tion 

S^. - j d'r 0;(r) Mr) , (5) 

where /x, v label the atomic orbitals, and Sfj_i, is the over- 
lap matrix between the orbital basis function and 
respectively. K^f^t is an additional fitting parame- 
ter with a value of 1.75 commonly used for molecules 
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and 2.3 for solidsii^iSi In the case of the planar 2D- 
graphene sheet a good match is achieved for a value of 
Keht = 2.8^ One important assumption within EHT 
is that the hopping matrix elements H^^^ , v fi de- 
pend linearly on the overlap matrix S^i, alone^ EHT- 
basis functions are usually Slater-type orbitals (STO), 
^nim = Nr^~^ e~^^ y/m(6,(p), which have the correct 
asymptotic form at large distance r (n, Z, and m denote 
respectively the principle, azimuthal and magnetic quan- 
tum numbers). The individual orbital wavcfunctions \q) 
are then approximated by a linear combination of STOs, 
with coefficients and exponents {ci,Ci} fitted for the in- 
dividual basis functions to match bandstructure data. 
Since the basis sets are directly fitted to experimental 
or theoretical data, the resulting AOs are more local- 
ized than the free atomic wavcfunctions, although they 
may still be regarded as long-ranged compared to the 
usual TB Wannier-like description; wc typically use a 
cut-off interatomic distance of i?c = 9 A for the inter- 
actions. The use of basis function that are not too lo- 
calized turns out to be a key issue for achieving a good 
transferabilityi2i Admittedly, the use of a direction in- 
dependent Keht function is a drastic assumption that 
can be relaxed by making the constant orbital depen- 
dent, but we make the simplifying assumption that the 
orientation dependence arises principally from the corre- 
sponding overlap functionals. 

The Hiickel prescription Eq. jsj can be generalized to 
connect different chemical sub-systems A, B. The prob- 
lem is that each sub-system has its own parametrization 
that yields its own vacuum level relative to which elec- 
tronic levels are calculated. For instance, the Fermi level 
of most metals is set to Ep = — 10 eV in the paper in^Si 
which means that each dispersion curve needs to be indi- 
vidually shifted to bring its Fermi level upto the experi- 
mental value, through the transformation H H + VcS 
for each sub-system. The correct alignment of the lev- 
els relative to each other becomes particularly impor- 
tant when studying compound systems such as metal- 
semiconductor junctions or molecular components at- 
tached on nanowires and nanotubes. We calculate the 
coupling matrix across such a hybrid junction between 
subsystems A and B using the Hiickel principle Eq.lO) 
through the interpolation scheme 

+ (XsHf, + y.s)] , (6) 

where VcA and Vcb are the shifts needed to align the vac- 
uum levels for sub-systems A and B, respectively. It is 
worth noting that this approach provides a simple inter- 
polation scheme that gives us the correct limiting val- 
ues of the Hamiltonian for the two individual subsys- 
tems. Further work, however, clearly needs to be done 
to calibrate this interpolation scheme to specific interfa- 
cial properties such as measured charge transfer doping, 
Schottky barrier heights, or perhaps first principles calcu- 
lations of intcrfacial dipoles or chemisorption properties 



such as workfunction modification at surfaces. 



IV. RESULTS: NANOTUBE 
BANDSTRUCTURES 

We begin by benchmarking our parameters with a two- 
dimensional graphene sheet. Figure |2] shows the corre- 
sponding dispersion relation along the M T K ^ 
M path within the 2D-Brillouin zone. The dashed line in 
Fig|21is the DFT-GGA calculation of the E~k dispersion 
calculated using the SIESTA code,S2i2^to which the EHT 
bandstructure is fitted (solid line) by adjusting the on- 
site Hamiltonian matrix-elements Hu^ the exponentials 
and the expansion coefficients of the Slater-orbital basis 
functions. As atomic like basis functions for each carbon 
atom we are using two basis sets: (i) sp^ and (ii) sp^d^ 
orbitals, each of which has been optimized to match the 
DFT-GGA bandstructure. The two parameter sets are 
given in Table ^ In our calculation all atoms within a 
cut off radius of 9 A from the two non-equivalent atoms 
within a unit cell are included. Furthermore, we set the 
Fermi level of graphene to Ep = 0.0 eV and A'cht is set 
to 2.8pS. 



C-graphite (EHT vs. siesta — - ) 




FIG. 2: 2D bandstructure of a graphene sheet along the 
M —> T ^ K ^ M path within the 1st Brillouin zone. 
The solid line corresponds to the EHT-bandstructure using 
parameters optimized to match the DFT-GGA bandstruc- 
ture (dashed line)i— *— The C-C bonding distance is set to 
ac~c ~ 1-44 A, and the cut-ofF radius for the neighbor atoms 
is i?c = 9 A. The Fermi level is at Ep = 0.0 eV. 
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C2 


C-sp 


2s 


-20.316 


2.037 


0.741 








2p 


-13.670 


1.777 


0.640 


3.249 


0.412 


C-spd 


2s 


-19.889 


2.025 


0.764 








2p 


-13.080 


1.624 


0.272 


2.177 


0.739 




3d 


-2.046 


1.194 


0.491 







TABLE I: EHT parameters for carbon fitted to tfie 2D- 
bandstructure of graphene calculated using DFT-GGA^2i— 
For botii parameter sets tiie parameter K^ht was set to 2.8. 



Carbon-nanotubes with chirality (n, to) are obtained 
by wrapping the two dimensional graphene sheet along 
specific circumference vectors C{n,m) = na.i -\- 
Note that this approach, being truly atomistic, goes be- 
yond the conventional zone-folding scheme. In the follow- 
ing, we will examine the transferability of the two EHT 
parameter sets (sp and spd) of graphene to bandstruc- 
tures for small diameter nanotubes. In our calculation, 
we assume that structural variations of the CNT affect 
its hopping elements Hij only through the overlap ma- 
trix Sij (cf. Eq. [SJ. This assumption means that the 
re-distribution of charge due to structural changes is dis- 
carded, so that the bandstructure is determined in a non 
self-consistent manner. For all tubes considered here we 
include the coupling of nearest neighbor unit cells which 
consist of two rings each in the case of armchair (n, n) 
tubes (translation vector Tg = 2.39 A), and four rings for 
zig-zag (n,0) CNTs (Tq = 4.32 A). 



A. Metallic armchair tubes 

Fig. 13 shows the bandstructure, transmission and 
density-of-states (DOS) for a (5, 5) armchair tube within 
a sp^-EHT model. The DOS shows typical features of a 
one-dimensional system, with a constant value within an 
energy interval of [—1.5, 1.5] eV, and van Hove singular- 
ities at the onsets of higher subbands. The transmission 
per spin within the interval [—1.5, 1.5] eV is 2 indicating 
that two bands can in principle contribute to transport. 
Including spin, one should then observe a maximum lin- 
ear response (zero bias) conductance of G = 2Go with 
Go = 2e^//i assuming no parasitic resistances posed by 
contact interfaces. Notably, curvature effects do not dis- 
rupt the bandstructure of armchair tubes, which stay 
metallic because the mirror symmetry is not broken upon 
wrapping the graphene sheet. 




Energy / [eV] 

FIG. 3: Bandstructure, transmission and DOS per spin for an 
armchair (5, 5) tube using sp^-orbitals. Due to their mirror 
symmetry the curvature of the tube does not break this sym- 
metry, so that the tube remains strictly metallic. The Fermi 
level is Ef ~ 0.0 eV and k is in units of [vr/To]. 



fcsT « 25 meVi^ at the Fermi energy. A simple 7r-orbital 
tight-binding model^S'^^'^i commonly employed in CNT 
transport simulations usually does not account for this 
effect, but such s-p hybridizations might be important to 
include when considering CNTs as possible candidates 
for interconnects. For tubes with diameters d < 1 nm, 
the structural deformation of the graphene sheet affects 
its electronic structure significantly enough that such an 
opening of a bandgap can be induced. The opening arises 
due to a reduction of the overlap between the nearest- 
neighbor TT-orbitals under deformation causing the Fermi 
wavevector kp to move away from the K-point within the 
1st Brillouin zoneM^ 

The more complex sp'^- and the sp^d^-EWT models we 
are using naturally account for these structural deforma- 
tions (Figs. 01 and [SJ through the structure-dependent 
overlap matrix S, cf. Eq.lQ, yielding for our bulk 
graphene parameters a gap that actually compares quite 
well quantitatively in the case of a sp'^-model with experi- 
mental scanning tunneling spectra (cf. Tablelll)l. whereas 
our sp^ d^ -model shows a poorer quantitative match for 
these class of tubes. 



Ultrasmall Diameter Tubes 



B. Curvature effects on non-armchair 
'semi-metallic' tubes 

Experimentally it is known that small diameter car- 
bon nanotubes that are normally expected to be metal- 
lic by the to — n rulcr^ with (9,0) and (12,0) chiral- 
ities for example, exhibit a curvature-induced gap than 



For ultrasmall nanotube diameters, hybridization ef- 
fects start to become dominant. Fig. |H1 shows that 
the (5, 0) zig-zag tube, semi-conducting in a zone-folding 
method, is predicted by EHT to become metallic for the 
case of a sp-orbital model, since the valence and conduc- 
tion bands cross at Ep ~ 0.0 eV. On the other hand, 
the zig-zag (6,0) tube moves the other way (Fig. [7J), 
from being metallic in a zone-folding method to semi- 
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Energy (eV) 
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FIG. 4: Bandstructure, DOS, and transmission per spin for 
a "metallic" zig-zag (12,0) tube using sp''-orbitals. The gap 
close to the Fermi level at Ep = 0.0 eV is about 50 meV. For 
comparison the experimental STS-dl/dV signal is shown at 
the bottom—. The bottom figure has been "Reprinted (ab- 
stracted/excerpted) with permission from M. Ouyang, J.L. 
Huang, C.L. Cheung, and CM. Lieber, "Energy Gaps in 
Metallic Single-Walled Carbon Nanotubes," Science, 292, 702 
(2001). Copyright 2001 AASS." 
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FIG. 5: Bandstructure, DOS, and transmission per spin for a 
"metallic" zig-zag (9, 0) tube using sp^-orbitals. The gap close 
to the Fermi level at Ep = 0.0 eV is about 80meV similar to 
STS-dl /dV measurements as shown at the bottom—. The 
bottom figure has been "Reprinted (abstracted/excerpted) 
with permission from M. Ouyang, J.L. Huang, C.L. Cheung, 
and CM. Lieber, "Energy Gaps in Metallic Single- Walled 
Carbon Nanotubes," Science, 292, 702 (2001). Copyright 
2001 AASS." 



conducting in EHT with a bandgap of « 0.12 cV due 
to strong hybridization effects. DFT-GGA calculations, 
however, show that a (6, 0)-tube remains metallic due to 
re-hybridizationi2Si^ Using in turn spd-orbitals for car- 
bon, our EHT-model agrees with the DFT-GGA results 
for these small diameter tubes. While the accuracy of 
DFT for semiconducting bandgaps is itself open to ques- 
tion, the contradictory result for the sp-orbital model 
could also imply that the inclusion of deformation effects 
through just the off-diagonal EHT Hamiltonian elements 
is no longer a valid assumption, and the onsite energies 
themselves need to be recalculated self-consistently to in- 
clude the effect of the electrostatics on the local atomic 
potentials. 

For large deformations a fully self-consistent calcula- 
tion of the bandstructure might be necessary to describe 
the electronic structure properly. If the tube has a large 
curvature, the respective charge distributions inside and 
outside the tube become asymmetrio^ leading to the for- 
mation of dipoles. Due to the charge re-distribution and 
the dipolar electric fields the individual bands are shifted 
such that the (6, 0) tube remains metallic in DFT-GGA. 
The processes of charge-redistribution, dipole-formation, 




Energy / [eV] 



FIG. 6: Bandstructure, DOS, and transmission per spin for a 
zig-zag (5, 0) tube using sp^-orbitals. Due to the crossing of 
the bands around Ep « 0.0 eV the tube becomes metallic. 
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{n, m) 


TB 


CNT-bands 


DFT 


Experiments 


EHT-sp 


EHT-spd 


(5,0) 




1.91 (7 = 2.5) 
2.29 (7 = 3.0) 


TTTTTT 




-0.05 


0.0 


(6,0) 


O.OSia 
f«0.3S 
0.179^ 


0.0 


-0.83^ 
O.O^i 


- 


0.12 


0.0 


(9,0) 


o.aitt 

0.075^ 

O.Oli^ 

^ 0.04^ 


0.0 


0.2 (GGA)ii 
0.17 (LDA)ffi 
0.12 (LDA)i*i 
0.17 (GW)S 


0.080 ± O.OOSii. 


0.075 


0.13 


(10,0) 


0.65^ 
0.8?^ 
0.85«i 


0.91 (7 = 2.5) 
1.09 (7 = 3.0) 


0.88 (GGA)iii 
0.8 {GGAp. 


0.83^ 


0.91 


0.95 


(12,0) 


0.0^ 


0.0 


0.08 (GGA)^ 
0.057 (LDA)^ 


0.042 ± 0.004ii 


0.045 


0.077 


(13,0) 




0.70 (7 = 2.5) 

fi C/i _ Q n-i 


0.73 (GGA)iii 


- 


0.71 


0.74 


(15,0) 


O.oai 


0.0 


0.14 {GGA)i^ 
0.038 (LDA)^ 


0.029 ± 0.004ii 


0.026 


0.05 


(16,0) 




0.56 (7 = 2.5) 
0.68 (7 = 3.0) 


0.61 (GGA)iii 


0.65±0.30ai 


0.59 


0.6 



TABLE II: Comparison of experimentally and theoretically determined bandgaps (in units [eV^]) using different theoretical 
methods: TB denotes orthogonal tight-binding, CNT-bands refers to a simple 7r-orbital description with one hopping parameter 
7, DFT is density-functional theory using different approximations for exchange-correlation potential, and EHT with sp- or 
spd-orbitals for carbons- 



studied tubes both with respect to STS experiments as 
well as other theoretical calculations. Our results based 
on Extcnded-Huckel Theory arc in good agreement with 
more sophisticated DFT approaches and experiments - 
in fact, local density approximation (LDA) suffers from 
the well-known bandgap problem that an EHT approach 
calibrated to graphene seems to bypass. This suggests 
that an EHT scheme, supplemented by self-consistent 
electrostatics, might be a good compromise between sim- 
ple TT-orbital tight-binding theories and computationally 
expensive DFT-mcthods, the practicality of the scheme 
being particularly important when modeling transport 
through large complex nanostructures. 



V. NANOTUBES AS CHEMICAL SENSORS 

While the previous examples test our EHT prescrip- 
tion for bandstructure calculations of bare CNTs, we 
now combine it with molecule chemistry. It has been 
suggested that the chemisorption of CO molecules could 
be enhanced by deforming the CNT, so that the re- 
gions of highest curvature have an enhanced chemical 
reactivity^ Based on a first principles calculation within 
GGA, da Silva et al. have shown that a semiconducting 
(8, 0) CNT can become metallic upon lateral deformation 
to an elliptical shape, but thereafter the semi-conducting 
state can be recovered by attaching a CO molecule at the 
spot of highest curvature iS4 This seems to indicate that 
highly deformed CNTs are possible candidates for sensing 
CO molecules by means of a bandgap variation and cor- 




-1 

Energy / [eV] 

FIG. 7: Bandstructure, DOS, and transmission per spin for 
a zig-zag (6, 0) tube using sp^-orbitals. The "metallic" tube 
becomes semi-conducting with a bandgap of ~ 0.1 eV, and 
shift of the Fermi-level from 0.0 eV to « 0.3 eV. 



and the floating of the bands, however, require a fully 
self-consistent bandstructure calculation that can be set 
up within the Poisson NEGF approach^i^ but has not 
yet been incorporated within our present computational 
scheme. Note, that this is not a shortcoming of the 
Extended Hiickel Theory itself, but arises instead from 
the importance of self-consistent electrostatic effects that 
have been discarded for simplicity in the bandstructure 
calculation. Tabic ^ compares the bandgaps for the 
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rcspndong metal-insulator transition. We use the DFT 
calculations by da Silva et al. as a qualitative benchmark 
to explore the accuracy of EHT for electronic structure 
in the presence of a periodic array of CO molecules, of 
which one unit cell is shown in Fig. |S| 

Our starting point is the bandstructure of an unde- 
formed semiconducting (8, 0) tube using spd-EHT pa- 
rameters optimized for 2D-graphenc with the Fcrmi-levcl 
set at Ep = —13 eVi2^ The transferability of the sp- as 
well as the spd-orbital EHT-parameters have been dis- 
cussed in the previous section. Setting our vacuum level 
as the zero energy reference, we shift the CNT Fermi en- 
ergy by 14 = -1-8.5 eV towards the experimental Fermi- 
level of 2D-graphene, i.e. Ep ~ —4.5 eV, employing 
the modified Hiickel prescription, cf. Eq. JHJ. Following 
Ref^ we deform the tube perpendicular to its axis so 
that the minor axis is 30% of the original tube radius of 
Rt = 6.3 A, cf. Fig.ia 



I 




FIG. 8: Unit cell of a semi-conducting (8, 0) CNT with one 
CO molecule at distance of d = 1.85 A from the CNT surface. 
The CO molecule has been placed above the center of the 
hexagon which is the most stable equilibrium position after 
relaxation—. The unit cell of the periodic structure consist of 
two (8, 0) CNT unitcells and one CO molecule, i.e. the cell 
contains 66 atoms. 

To ensure that our strongly deformed structure really 
corresponds to the local minimum of the total energy, 
we optimize the deformed tube unit cell by means of 
Gaussian 03^ using the spin unrestricted LDA approach 
within the Vosko-Wilk and Nusair (SVWN) approxima- 
tion. During the optimization we froze two rows atoms 
along the opposite extremes of the minor axis, while edge 
effects were eliminated from the optimized structure by 
employing periodic boundary conditions within Gaussian 
03 for structure optimization, with a translation vector 
of 4.32 A along the tube axis. We note that compared to 
da Silva et. aU^ our system being optimized consists of 
only one unit cell of a (8,0) tube instead of two. 

Since the EHT levels are correct upto an overall shift, 
we start by aligning the levels of the CO molecule and the 
CNT. The highest occupied molecular orbital (HOMO) 
of CO, calculated within Gaussian 03 using the Becke 
3 parameter exchange with Perdew-Wang 91 correlation 
(B3PW91) with 6-31g(d) basis sets is Eh = -10.14 eV 
relative to vacuum. The corresponding value in our EHT 
parametrization is Eh = —14.09 eV, so that we need an 
overall shift Vc^co = +3.95 cV. After this shift is exe- 



cuted, the CO molecule is attached to the tube at the 
point where its curvature is highest. According to RefiS4 
the most favorable location to place the CO molecule 
is above the center of a hexagon, as shown in Fig. |S1 
Since we are only interested in the effect of one single 
CO molecule on the tube dispersion, we need to avoid 
any overlap between neighboring CO molecules. This is 
accomplished by using two CNT (8, 0) unit cells and at- 
taching only one CO molecule. The effective periodicity 
of the system is then 8.64 A, large enough that the neigh- 
boring CO molecular basis functions do not overlap. 

Figure |H1 shows the dispersion, and Fig. ^| the re- 
spective DOS and transmission per spin for the com- 
bined CNT-molecular system. The original bandgap 
Ea = 1.1 cV of the undeformcd tube (left) shrinks down 
to Eq ~ 30 meV upon the 30% lateral deformation, so 
that the tube becomes effectively metallic at room tem- 
perature (center). Attachment of the CO molecule on the 
deformed CNT (right) makes the tube semiconducting 
once again, with a bandgap of Eq « 100 meV, much less 
than that of the undeformed tube, but noticably larger 
than thermal energies at room temperature. Our results 
agree qualitatively with da Silva et. al^^i^ even though 
there are quantitative differences: the initial bandgap 
of the undeformed tube is 0.66 eV and becomes com- 
pletely closed after lateral deformation, so that the tube 
becomes truly metallic. The recovered bandgap upon 
CO-attachment is with « 200 meV, of similar order as in 
our case with 100 meV. The largest source of disagree- 
ment is the bandgap of the undeformcd CNT; given that 
DFT bandgaps (see Table |nj are often questionable as 
well, a quantitative resolution of this discrepancy may 
need experimental attention. 



undeformed deformed deformed + CO 




FIG. 9: Bandstructures of a zig-zag (8, 0) CNT using sp^d^- 
orbitals. Left: undeformed tube with a bandgap of Eq ~ 1.1 
eV, middle: deformed tube with a small gap of Eq ~ 25 
meV, and right: deformed tube with attached CO molecule 
at distance 1.85 A. The gap here is Eg ~ 100 meV. The Fermi 
energy is indicated by the dashed line. 

Note that the present calculations arc all non-self- 
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Energy / [eV] 

FIG. 10: DOS and transmission per spin for the semi- 
conducting (8, 0) CNT: left: undefomed, middle: deformed, 
and right: deformed with attached CO-molecule. 



consistent, ignoring effects due to the rearrangement of 
charges during tube deformation, as well as those aris- 
ing from charge transfer between the molecular species 
and the nanotube, while a proper self-consistent calcu- 
lation as discussed in the last section, is needed to do 
quantitative justice to this problem. Nonetheless, the 
EHT parameters seem to be quite transferable between 
bandstructure as well as surface chemistry, in particu- 
lar for strongly deformed structures without the need 
for any reparametrization. This makes EHT a good 
compromise between accuracy and practicality. In our 
follow-up paper, we will demonstrate the applicability of 
EHT to modeling silicon, including its transferability be- 
tween bulk and multiple surface bandstructurcs of recon- 
structed silicon surfaces and also for nanowircs. Apply- 
ing this aproach has allowed us to quantitatively explain 
and in some cases even predict interesting experimental 
results involving molecular conductors on siliconi^ 



VI. FUTURE WORK 

For transport calculations, we often need a minimal 
model that can do justice to bandstructure, electrostatic 



as well as bonding chemistry. This becomes particularly 
important if one wants to deal with strongly deformed 
structures, interfaces or combinations thereof including 
molecules. We have shown that Extended Hiickel The- 
ory provides a good practical compromise to capture var- 
ious aspects of bandstructure and chemistry. The two 
attributes that make this semi-empirical approach espe- 
cially useful are the presence of explicit basis-sets and 
non-orthogonality. To make this chemists' approach to 
bandstructure of further use, it may be preferable to work 
with a self-consistent complete neglect of differential over- 
lap (CNDO) approach to bring in differential Coulomb 
costs into the picture. 

For many nanoscalc structures such as nanotubcs 
and nanowircs perhaps even interfacing with smaller 
molecules, we believe that a semi-empirical approach 
combining bandstructure and chemistry is essential, 
given that typical tight-binding theories are not trans- 
ferable beyond small deformations while DFT theories 
are computationally quite prohibitive. The latter be- 
comes even more difficult to implement when we want to 
move from equilibrium electronic structure to more com- 
plicated nonequilibrium transport problems. It is gener- 
ally acknowledged that much of the conducting proper- 
tics of nanostructures are dominated by their interfaces 
and contacts. It is in this complicated domain that we 
believe the real strength of a non-orthogonal theory with 
explicit basis sets such as EHT or CNDO is likely to 
manifest itself. 
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